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Annotation 


The paper considers the problem of energy and multicomponent 
ambipolar diffusion of plasma in the lower ionosphere of a planet 
with a weak magnetic field. Energy and diffusion equations are de- 
rived for three-temperature plasma in which heat and mass transfer 
processes are taking place jointly with chemical reactions in a 
form which is convenient for calculating models of the composition 
of a multicomponent ionosphere. Unlike in some early formulations 
of the problem [1-4], it is shown that the coupling effect of the 
electric polarization field on the motion of charged particles in 
the ionosphere is such that the diffusion rates of individual ion 
components are determined in terms of number density gradients of 
all ion components of the mixture. The proposed approach general- 
izes the known description of ambipolar diffusion derived earlier 
as part of the theory of a traditional ternary mixture (e,i,n) [5-9] 
to the case of an atmosphere composed of many charged components. 
Energy interactions among components are analyzed, in particular, 
thermal energy losses or electrons (ions, neutral particles) due 
to elastic and inelastic collisions, which allowed to obtain energy 
equations for each component of the system. 


MULTICOMPONENT DIFFUSION AND ENERGY C1L rtACTERISTICS OF PARTIALLY 
IONIZED PLASMA IN THE IONOSPHERE OF A PLANET 

M. Ya. Marov and A. V. Kolesnichenko 

Introduction 


The transfer of charged particles by means of diffusion, which 
along with the photochemical process controls the vertical distri- 
bution of ionosphere components is an Important factor in all diverse 

i 

aeronomy processes taking place in the ionosphere of a planet. The 
diffusion of charged particles in weakly charged multicomponent 
plasma of the lower ionosphere differs from diffusion of small com- 
ponents in the conditions of a neutral atmosphere. This is related, 
in the first place, to the necessity of taking into account electro- 
static forces not only during the collision of charged particles, 
but also during the collision of charged and neutral particles, 
when an electric dipole moment is induced in the latter, and in the 
second place, to the necessity of taking into account the electric 
field effect of a space charge formed during the fast diffusion 
movement of electrons compared with the motion of ions, impeding the 
relative diffusion of charged particles. 

In the case of three-component plasma (e,i,n), the electron- 
ion plasma undergoes joint (ambipolar) diffusion under the effect 
of an electric polarization field 2 with a common diffusion coeffi- 
cient and common velocity Problems of plasma diffusion in 

the ionosphere assuming a ternary mixture are discussed in detail 
in a number of studies (see for example [5,8], and also the joint 
authorship study [9]). 

For weakly Ionized ionosphere plasma, consisting of electrons, 
various kinds of ions and neutral partiules, the diffusion rates of 
individual charged components are not equal to each other and are 
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determined from the condition which states that the longitudinal 
electric current is approximately aero. Few studies discuss the 
ambipolar diffusion problem in the general case of a multicomponent 
ionosphere (see for instance, [1-4, 9-11]), while several studies 
along these lines do not use equivalent expressions for diffusion 
rates of Individual ion components of the gas mixture (for Instance 
[1] and [4], 

The purpose of this study is to derive energy equations for 
multi-temperature plasma of the lower ionosphere and also general 
diffusion equations for determining diffusion flows for types of 
ions isolated in the mixture of atmospheric ions using the concept /5 
of effective ambipolar diffusion coefficients. Another objective 
of the study is to eliminate a number of inaccuracies in the de- 
scription of ambipolar diffusion of a multicomponent conductive 
medium which occurred in some early formulations of the problem 
[ 1 , 3 , 4 ]. 

1. Stefan-Maxwell Relations for Multicomponent Diffusion 


For a quantitative analysis of multicomponent diffusion pro- 
cesses in ionosphere plasma, we will consider the Stefan-Maxwell 
relations derived in classical kinetic theory of multicomponent 
gaseous mixtures of monoatomlc gases of moderate density as part of 
a first approximation for multicomponent diffusion coefficients and 
a second approximation for thermal diffusion coefficients [12] 




r® $( 2 h 
1 S K 5, J* 
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Here m , p , n is the molecular weight, mass density and number 
s s s 

density of the s-th component, respectively, V , x -n /n is the 

s s s 

diffusion velocity and molar concentration of the s-th component, 

k, p, T is the total number density, mass density and temperature 

T 

of the gas mixture and , D gk are respectively the thermal diffusion 
coefficients and binary diffusion coefficients for all pairs of 


2 



components. The vectors of diffusion forces a s are defined as 



where c s =p g /p Is the mass concentration of the s-th component, p 
Is the pressure of the gas mixture and F is the force acting on the 
s-th component. The order of approximation 5=1,2,... with which the 
transfer coefficents are determined In (1.1) corresponds to the 
number of first terms in a series expansion of the coefficients of 
perturbed distribution functions of components In terms of Sonine's 
polynomials . 

Until recently, relations (1.1) for diffusion rates were only 
derived in first approximation in the kine.ic theory of gases. Usu- /6 
ally such approximation is considered to be adequate, however in 
ionized gases, the admissible error may turn out to be considerable. 

The necessity of taking into consideration higher order approxima- 
tions in calculations of transfer coefficients in an investigation 
of flows of ionized mixtures of gases in the ionosphere (5=3 and 
higher [12,13]), when not only the temperature and pressure but also 
the elementary chemical composition [14] of the flow undergoes 
changes, makes the Stefan-Maxwell equations in form (3.1) inappli- 
cable in aeronomy in the general case. 

However, study [15] gave very recently a derivation of relations 
(1.1) from the kinetic theory of gases in any approximation for trans- 
fer coefficients using the Chapman-Enskog method, and study [16] de- 
rived relations (l.l) for the case of imperfect (perfect) mixtures 
of gases by methods of thermodynamics of irreversible processes, 
while demonstrating the symmetry of the matrix of resistance to dif- 
fusion, in complete agreement with the latest results of the kinetic 
theory of gases. 

Let us make some transformations. Henceforth we will consider 
as external forces electric, magnetic and gravitational forces 
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where $ Is acceleration due to gravity, e is the electric charge of 

s 

the s-th component (for neutral components we will assume e HO), 

. s 

E is the electric field Intensity vector, B Is the magnetic Induction 
vector, which In the general case Includes the unperturbed magnetic 
field of the planet and a small induced field. Using the relation 
P a =x s p for the partial pressure of the s-th component, we have 
Vx ■■Vp /p-x Vlnp, from which we obtain 

o S o 

(1 - 1 ° 

where e n v is the density of the total conduction current in 

5 S 5 5 

the plasma. 

Expression (1.4) takes into account the quasineutrality of 
ionosphere plasma 

H 

Ze 4 ti s *0 ( 1 . 5 ) 

Using formula (1.4) and also the full equation of motion (without 

the viscous term) for the continuum which models the mixture of gases 

as an entity, we rewrite the Stefan-Maxwell equations (1.1) In the 

form of equations of motion for individual charged components of /7 
1 

the Ionosphere , 

3.5J (1.6) 

In the above, v g j=KT/nD s j are the coefficients of friction which are 
either determined experimentally or from a detailed investigation of 
the dynamics of particle collisions using the kinetic theory of gases, 
K is Boltzmann’s constant, N’ is the number of charged components 


1 The viscosity of the system is almost completely determined by the 
neutral component of the ionosphere [ 7 ] . 




of the gas mixture, and a is the thermal diffusion factor defined 

s 

by 



which constitutes a measure of the relative importance of thermal 
diffusion and ordinary diffusion. In the general case, the thermal 
diffusion factor Is a complicated function of temperature, concen- 
tration, and molecular weights depending parametrically on the laws 
obeyed by intermolecular forces [12], According to available esti- 
mates, thermal diffusion for principal components of the ionosphere 
is a negligibly small quantity due to the small difference in mole- 
cular weights, which essentially should only be taken into account 
in light gases (hydrogen and helium [17]). 

Equation (1.6) did not take into account the possible aniso- 
tropy of coefficients of friction in the magnetic field, which is 
related to not having taken into account the effect of the magnetic 
field in Stefan-Maxwell equation (1.1). At the same time. It is 

clear that In strongly magnetized plasma (ft ft = ), in 

6 61 6 

particular, such as the ionosphere in the P region of the 'Earth’s 
ionosphere, the collision term must be anisotropic. In simple cases 
(for example for completely ionized two-component plasma), in ion- 
osphere conditions, this effect can be estimated by the order of mag- 
nitude [18]. We also note that although the presented derivation of 
equations of motion (1.6) is only valid in a single-temperature ap- 
proximation, the form of the equations remains unchanged also in the 
case of multi-temperature, plasma [19]. 

t 

Directly in the ionosphere of the planet, neutral molecules are 
much more numerous than charged particles: n n >> n g , n^. Therefore 
in describing physical processes in the lower ionosphere, one can /8 
often restrict the analysis to the theory of a weakly ionized gas. 

In this case, v/e can assume that the Hydrodynamic velocity coincides 
with the velocity of the neutral gas , all components of which 


5 



have the same velocity V (assuming gravitational-diffusion 
equilibrium in the thermosphere) . Then the general eouation 
of motion for the plasma component of the cas mixture, 
obtained by summing equations (1.6) over all charged gas com- 
ponents, assumes the form 



( 1 . 9 ) 


Here \i Dn is the total coefficient of friction of the s-th 
charged component with the neutral component of the ionosphere ; 
n n is the total number density of neutral components of the 
gaseous mixture, p l , n L , P L are respectively the mass density, 
number density and pressure of the plasma component of 
the ionosphere. The upper index "n" in a sum denotes 
summation over neutral components of the system. The general 
equation of motion for the neutral component of the 
ionosphere 





( 1 . 10 ) 

( 1 . 11 ) 
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follows from equation (1.8), and the full equation of 
motion for the continuum modeling the gaseous mixture as 
a whole entity 


where II is a tensor of viscous stresses of the gaseous 
mixture, n n is the viscosity coefficient and p n , p n 
are the mass density and pressure of the neutral component 
of the ionosphere. 

2. Collision Effects 


Detailed calculations of various transfer processes 
in the ionosphere of a planet require a rigorous treatment 
of the problem of collisions among various particles. 

The collision term 

in equation (1.6) represents the velocity with which a 
unit of volume of the k-component loses momentum as a 
result of collisions of its particles with other components. 
Introducing the mean collision frequency S Jcs = n s v i < s rn— > 
we express the collision term in the form /9 

cA\ rts = n, ^ ( V* - » a A tn KS - vy . _ ( 2 * 1 ) 


7 


In the above m lcs =m j <: m s / C m j c +m s ) I s the reduced maos, 0kQ" n 6 Sks^ks is 
the r.ean collision frequency with momentum transfer between particles 
of type k and s, g kQ = [8KT R /Hm ks ] 1 ' /2 is the mean relative speed of 
particles with Maxwell velocity distributions, T R *=m kc ( T k + J s is the 

m k % 

reduced temperature and Q kg is the mean collision cross-section with 
momentum transfer. We shall present some estimates: The scatter- 

ing cress-section Q ks is determined from an experiment or by means 
of gas kinetics calculations [12,18,19], when the law governing the 

interaction between particles is known sufficiently well. Kinetic 
* 

theory gives the following expression for Q k£l in the case of collision 
of lens with a neutral gas [19] 




( 2 . 2 ) 


where Q in (g) Is the momentum transfer c.ross-sectlon between particles 
of type i and n, having relative thermal speed g and equal macro- 
scopic velocities of motion . A more exact expression for Q ln 

taking into account differences in the velocities of motion of com- 
ponents i and n is presented in study [20]. Rigorous analysis shows 
that as a charged particle approaches a neutral particle, an electric 
dipole moment is induced in the latter and the collision cross-section 
between the particles is determined to a considerable degree by the 
resultant electrostatic forces. For ions with neutral particles. 

Banks [21,22] proposes the following semi-empirical expression 


Qin-W -10'k 


*K(* + £JI‘ •« S 


(2.3) 


_ 2 1 | "< 

w her* the atomic polarisation a (=10 a' cm ) for the principal 
neutral components of Earth's ionosphere has the following values: 
0*^*1. 60, a» N =0.82, a' Q =0.89, a' M =0.67, ct' He =0.21, a' N =1.13 [23,24]. 

It is assumed that the ions and the gas obey a Maxwell distribution. 

In this case, the collision frequency is Independent of the temperature 

and is determined from relation [21,2] ORIGINAL PAGE IS 

OF POOR QUALITY 


8 


K< S ‘ 10 '\(</i “ a )'1 sec" 1 <2. H) 

nj f\ 

where [i and A Is the mass of the particle In atomic units. /IQ 


A,A„ 


The resonance charge transport reaction plays an Important part 

in the diffusion of ionB of some type in the Bame type of gas. The 

frequency of collisions for this process can be approximately repre- 

1 

santed in the form 

( Tc + T k ) n t CCK - 1 (2.5) 


2 1 /? 

/3ec*deg. ' ) for the prin- 


v/here the constants 3? ( e 10 _12 s! cm 

m m 

clpal components of the Earth's ionosphere have the following values: 
6 O'tcf 1 - 6 - s'oi.o/ 1 - 1 ’ 6 HtH =1-0 > 6 HetHe' 3-0 * S 'n 2 ,N 2 =2-1 ’ 


C293. 


While considering the collision frequency of electrons with a 

heavy particle, we will use the condition m <<rn ; then m - m , 

e n en e 3 

- T g , and we obtain for the mean collision frequency 



(2.6) 


For the principal neutral components of the Earth's ionosphere, the 
mean collision cross-sections with momentum transfer have the values: 
Vn = (3.76-4.54*10"** T e ) T e 1/2 -10“ 17 , cm 2 , Q q)(] = 2.93*10" l6 "(H-3.6* 

10- 1? .T e 1/2 ), cm 2 ; Q e , Q " (4.53-1.3) 10" 16 , cm 2 ; Q q ^ (7 2 . 9-9 . 93 • 107 3 Te, 

cm 2 ; Q e He = (7.46-0.8) 107 l6 cm 2 [21,22], 

The total collision frequency of an electron in multicomponent 


2 In all formulas applied in practice, the temperature is expressed in 
degrees Kelvin and all other magnitudes in CGS units. 
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plasma is equal to the sum of collision frequencies with different 
components 0 = £_ 0_„ . 

0 060 

The elastic collision dynamics of two charged particles are 
described by a well known potential, which allows to derive an ex- 
pression for the mean momentum transfer cross-section [18], 



e /\ct 




cm 2 (2.7) 


Where t is the dielectric permeability coefficient of the plasma and 
A e i= 1. 24 • lO^v^’^/n^ is the Coulomb logorlthm (ni). The correspond- 
ing mean electron collision frequency with momentum transfer is 


%r n i%iQ« n i 


fr* Rz/ fn* c a fytt'itGiAr,: 

3( XT,/ ’imjMv* 


sec 


-1 


( 2 . 8 ) 


For ion-ion collisions, according to [21,22], the cross-section /II 
can be expressed in the form 


Q il =^,^’fO' 6 en.AL}/T^ . , cm 2 (2.9) 

p 

where the parameter A., is defined as A. ,= 2E R n /e , E is the mean 

1J ±J U 2 * _ 

energy of the relative motion of two ions and R n = 4 re /Kc (ni/Tl + 

— ] /? ^ 
nj/Tj) is the Debye radius. It is assumed that both gases have 

a Maxwell velocity distribution. Then the corresponding frequency 

of ion-ion collisions is 

&j Qir 8 ’ 1 * 4 . (2.10) 


3, Amblpolar Diffusion of Ionosphere Plasma in the Absence of 
Ex ternal Electromagnetic Forces 

We will now consider diffusion processes in a multicomponent 
ionized gas in the ionosphere of a planet with a weak magnetic field 
§ £ 0, on the assumption that all different types of particles par- 
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ticlpating in this motion have a different temperature. Because of 
the smaller mass, electrons undergo faster diffusion than Ions, how- 
ever the space charge formed In the process gives rise to an electric • 
field £ which tends to slow down their motion. The electric field 
acts in the opposite direction on Ions and accelerates them. Ignoring 
the magnetic field Induced by the separation of charges, the Stefan- 
Maxwell relation (1.6) for a charged ionosphere component s can be 
expressed as follows 


Kl 7\ 3 - e, n s E ^ ^ £ J, ( s« U ... ») 

frit -> n, / * \ n 


(3.1) 

(3.2) 


Here J = n (£ -^) is the diffusion flow (flow of the number of par- 

o S X 

tides of type s with respect to the mean mass velocity V‘ S a v£ * 

t «(j/g is the unit direction vector 5 and H =KT /m n g is the local 

altitude scale of the s-th component. Introducing the diffusion 

movement coefficients C and the diffusion coefficients D of charged 

s s 

particles of type s 






k 


(3-3) 


relations (3.1) can be written in the form 






(3.1*) 


/12 


Relation (3.3) between B and C_ is the Einstein relation, known 

s s 

earlier for a ternary mixture of neutral components, ions and electrons. 


Let ns write down the expression for the electric current den- 


sity 




Taking into account (1.5) and (3.1), we have 


. TT t— i 




(3.4) 


ll 


where the conductivity of the plasma (longitudinal conductivity in 
the presence of a magnetic field) is determined from the expression 

N U 




(3.5) 


Since a minute vertical separation of electrons with respect to 
positive ions leads to the formation of a strong electric polarization 
field in the ionosphere impeding the relative diffusion of charged 
particles, ultimately a stationary state is established in ionosphere 
plasma, in the presence of which ions and electrons move in such a 
way that a space charge does not arise. In this case, only small 
deviations from equilibrium can occur and any vertical electric 
current must be negligible. Setting in expression (3.4) ^-0, we 

shall find the electric polarization field E h under the action of which 
ambipolar diffusion is taking place in the ionosphere. Next, elim- 
inating this field from relations (3-D and omitting the subscript 
h in the vertical component of vector quantities, we write the Stefan- 
Maxwell relations in the form 


a t) N .< 

KTi A, - «T i \ *&C(nJr ".*! «■»...». '•) (3.6) 


So-called ambipolar coefficients of friction ti.- ' 

were introduced here. In deriving formula (3.6), for symmetry of the 
final relations, the zero term e< j «2.2« —O' (by virtue of the 
quasilinearity of (1.5) and relations (3*3)) was added to the right 
member of (3.1). 


Relation (3.6) can be simplified. For atmospheres of planets, 
electron mobility exceeds many times ion mobility |C e |>>C 1 [25]. /13 

Consequently, the following approximate relations are valid 


K<?. 




J * 


WS * A = [ W A. T. i A.] - £ rfVW 
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and equation (3.6) for vertical diffusion of type 1 Ions can be ex- 
pressed as follows 


A, + £ KliA e -J «'■'( nj'-njj, (i-u...*!/ (3 7) 

V S’~ >*. * 'L rl e i n * K./& a »i, . 


In the case when electrons and Ions have the same temperature 
(T e “Ti=T), It is convenient to rewrite relation (3-7) In terms of 
amblpolar binary diffusion coefficients 




'k« 


KT 


1 G W 

=1.5;." 2 




cX>«* -1 + 


(3.8) 


in the form 



( t * ^ 2, . . , ^ 


(3.9) 


In investigating the composition of the atmosphere of a planet, re- 
lations (3.9) together with differential continuity equations for the 
concentrations n^ of ion components of the atmosphere 


dL 

c a 


( ft* /s) J, lif- Ki II ", a ] , 


( 3 . 10 ) 


and the plasma quasilinearity condition 

(3.11) 

allow to calculate the vertical distribution of charged components in 
the lower ionosphere of the planet at a given temperature profile. 
Here n ir , Y ir are the stoichiometric coefficients in the r-th (r=l,2, 
. . . ,R) aeronomie reaction, *]*... “k'* &*W + ?a*ftK.. ; ft is the 

number of reactions, Kr, Kr’ are the constant velocities of the 
direct and reverse reaction respectively. 
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Equations (3.9) generalize to the case of diffusion of various 
types of Ions and electrons in a mixture composed of many components 
the well known diffusion equations introduced by Kolegrov et al [26] 
for determining the densities of individual components used in many /l^ 
studies during calculations of models of the composition of a neutral 
atmosphere. 


Ignoring for simplicity the inertial term in expression (3.2) 
for A , and restricting ourselves to a consideration of stationary 

O 

movement of a neutral atmosphere, we write relation (3.7) in terms 
of the gradients of principal hydrodynamic quantities in the fol- 
lowing final form 


vTt Vr t Vk ^ c+ vn ~ 

23 vS(nj K -iijj/Kii t (T € +Ti) # 


I**. 


1+oio 


(3.12) 


where (T i +T e )/m 1 g is the altitude of the homogeneous atmosphere 

of the i-th ion component. The influence of the thermal diffusion 
effect on the altitude distribution of ion components of the atmos- 
phere is only considerable for minor components [27]. Equations 
(3.12) for diffusion flows of ion components are fundamental in a 
descritpion of ambipolar diffusion and generalize the results obtained 
earlier [5-10] for three-component plasma to the general case of 
weakly ionized multicomponent plasma in the lower ionosphere. 

Finally, assuming that collisions among ions and neutral par- 
ticles play a greater part than collisions among charged particles 
(which is valid for the lower ionosphere in the E and F regions ) i 
and also taking into account the fact that in the thermosphere the 
neutral components are in a state of gravitational-diffusion equili- 
brium (Jj_ s 0), we obtain from equations (3-7) an expression for the 
diffusion rate of the j-th principal plasma component 


^ (3.i3) 
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where the ambipolar diffusion coefficient is defined by the 

relation 


ay. rL _ K(Tc+m& , H) 

J KfVir TjD^TcDj " 


The index n in the sum denotes summation over neutral components of 
the system. Thus, on the strength of relation (3.11), the drift vel- 
ocity of each Ion component depends on the gradients of all ion com- 
ponents of the ionosphere. 


In the special case of three-component plasma (n e =n i =n), form- 
ulas (3.13) and (3.14) become the well known relations for the coeffi- 
cient and ambipolar diffusion rate of plasma in a coordinate system /15 
moving with the velocity of the entire substance 




4. Ambipolar Diffusion in Proper Magnetic Field of Planet 

In discussing the problem of diffusion in the Ionosphere of a 
planet with its own magnetic field £, we will limit ourselves to the 
case in which It suffices to consider the Interaction of charged 
ionosphere components only with neutral particles (friction with 
neutral gas) which are in a state of diffusion equilibrium, without 
taking into account collisions of charged particles among themselves. 
In this case, the equation of motion (1.6) for each type of charged 
particle will be represented In the form 




K u 


OftGHVAL ,-„ uE , s 
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Denoting by t jj and the components of an arbitrary vector 
along and across the line of force of the magnetic field IT / 

h.=B/B ) , we write the longitudinal and transverse com- 
ponent of equation (4.1) in the form 



/ 


(4.2) 


(4.3) 


where ft . = e B/m c is the generalized gyromagnetlc frequency of the 

s •K 3 s K ± 1 d- 

s-th component, and E r = E + ^ V x B Is the electric field in a 
coordinate system moving at the velocity of the neutral gas in the 
ionosphere . 

•J-'.ing formulas (4.2), let us write down the expressions for the 
density of the flow along and across the magnetic field 


i 

\K. 


T' T* 


(4.4) 


where the transverse and Hall conductivity is determined from the /16 
expressions 


~ Jj e s rl Aj. , 4s = Jj e i ^ ^sa 


(4.5) 


In Earth's lower ionosphere (F-reglon) B s 0.3 gauss, 3., ..Msec - ^ 

3 35 sec - , 3 .^1500 sec", m 0 s 9.108*10~ g, e a 1.602*10“ e.m.u. 

an si e 

and the following estimates 

«i, KMM (li S) 
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are valid in a considerable range of altitudes. Using these estim- 
ates, we obtain the following approximate relations 




(4.7) 


Eliminating the longitudinal electric field 2^ from relations 
(4.2) and (4.4) and ignoring terms which are small compared to unity, 
we obtain the following expression for the longltudina.1 diffusion 
flow of ions 



The term involving J., can be omitted. Indeed, C n /C n>10"? 

— 12 ' g S S G S 

J a *10“ e.m.u. [ 28 ], which together with n^lO cm, gives iJ^/en |* 

P 1 I 

10 cm/sec. Since (V^ j in the ionosphere is probably of the order 
of magnitude 1 m/sec, the term involving 2 in relation (4.8) is 
small, so that we can assume 



(4.9) 


In complete accordance with formula (3.13). 


To obtain the transverse component of the diffusion Ion flow 
In ionosphere conditions (4.7) , we proceed as follows: Using the 

full equation of motion of the plasma component of gas mixture (1.8) 
in the form 



(4.10) 
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we write the last term in the right member of expression (4.4) for 
as follows 



(4.11) M 


Then from (4.4) we obtain for the transverse component of the 
electric field 






(4.12) 


Next, eliminating SJ_ from relations (4.2) and (4.12) and omit- 
ting terms which are small by virtue of conditions in the ionosphere 
[4.6], we obtain for the transverse diffusion ion flow the following 
final expression 






(4.13) 


in which the last term in brackets is eliminated using the full 
equation of motion of the neutral gas mixture component (1.10) 


h' n 


l tx =3 <» ■ » 


Formula (4.13) generalizes to the general case of weakly ion- 
ized multicomponent plasma the results obtained earlier for three- 
component plasma [6,7]. Ignoring for simplicity the inertia of 
particles and also the effect of thermal diffusion in equation (4.l4), 
we obtain from (4.13) in the special case of a ternary mixture 
(C^n^/fi^l/e ) the following standard expression for the transverse 
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dil fusion rate of electrons (Ions) [6]. 



£(^sTr«pPe>E\ 


(4.15) 


If r s denotes the rate at which the s~th component of the ion- 
osphere is formed as a result of all aeronomic reactions, a good 
approximation of the continuity equation for the s-th ionic compon- 
ent i.i the following equation 


i. 

c a 




K(r c ^) 



(4.16) 


where J g is determined from formula (4.13), which allows to cal- 
culate the distribution of charged components in the lower ionosphere 
of a planet with a magnetic field when the temperature profiles in /18 
the ionosphere and also the dynamics of the neutral gas are given. 

The degree of complexity of the right member of equation (4.16) de- 
pends on which reactions and components are acknowledged to be neces- 
sary for a description of properties of the ionosphere at the given 
altitude. Equations (3.12) and (4.16) derived in this study are 
promising for a description of ambipolar multicomponent diffusion of 
Individual types of ions in plasma in the lower ionosphere. Numer- 
ical modeling of the vertical structure of the ionosphere based on 
these equations is of considerable interest. 

5 . Deviations f r om Local Thermodyna m ic Equilibrium in the Ionosphere 
Of a Planet ~ ’ 


A characteristic peculiarity of the ionosphere of a planet is 
the absence of thermal equilibrium, which is reflected, for example, 
in a difference in kinetic temperatures of electrons, ions and 
neutral particles. The basic reason for the absence of local 
thermodynamic equilibrium in the ionosphere is the circumstance that 
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In the general case, the energy of solar photons or charged space 
particles is greater than the energy necessary for ionization. A 
great portion of this energy is carried away by photoelectrons. 

Their initial energy spectrum 

(s.d 

is a complicated function of an incoming photon hv, the ionization 
potential of an atom subjected to the effect of radiation 

and also the degree of excitation (Dj is the excitation energy of 
the j-th component to the a level) of the positive ion that is 
being formed. The energy of fast photoelectrons goes partially into 
heating the electron gas thermallzed by elastic collisions, into 
excitation of electronic vibrational and rotational energy levels of 
molecules and the fine structure 0(^p), and partially into formation 
of excited ions and secondary electrons, which in turn are capable 
of exciting electronic and vibrational energy levels of atmospheric 
particles and elastic interactions with thermallzed electrons. The 
heating of electron gas in the ionosphere during collisions with 
photoelectrons and secondary electrons leads rapidly to a Maxwell 
distribution of velocities characterized by the electron temperature 
T (sometimes "tails" of high-energy superthermal electrons may /19 

exist in the upper ionosphere and protosphere, since collision!? at 
high altitudes are rare, so that during their existence the electrons 
are not able to undergo thermallzation) . The heated electron iono- 
sphere gas which is generated is cooled on account of elastic and 
inelastic collisions with neutral particles having temperature T^ 
and also during Coulomb interaction with ions in the surrounding 
medium having temperature T^. In principle, different ion components 
may have different temperatures (for example, in the Earth's ion- 
osphere, the temperatures of light and heavy ions may differ by 100°K 
[30,31]); however in the general case, the latter is impeded by the 
Coulomb interaction between them. Heated ionosphere ions are then 
cooled in elastic collisions with neutral particles. 
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We will present a. quantitative description of the process of 
equalisation of the temperatures of Ions and electrons In the Ion- 
osphere, while limiting ourselves for simplicity to the case of 
two-component plasma. We denote by T the time at which the Maxwell 
distribution of electrons Is established (thermalization time), which 
occurs as a result of Internal Interactions among the elections pro- 
per. The analogous quantity for Ions will be denoted by T,.. Pin- 

£ XI 

ally will denote the energy exchange time betv/een electrons and 
ions (time in which local thermodynamic equilibrium Is established 
within the system) . 

Rigorous theory [32] gives the following expressions for the 

times T , T 

ee el 


WareX&iA < ■ (5-2) 

Here In A is the Coulomb logorithm (which varies from 10 to 20 for 
typical conditions in the ionospheres of planets). If both elec- 
trons and ions have approximately a Maxwell distribution with re- 
spect to velocities with temperatures T e and T^, the process of 
equalization of the electron and ion temperatures is described by 
the relaxation equation [33] (see also equation (6.11)) 


dT* Te-H 

~~ Tft J 


= JL m m \ ]• Tc 1 %. 


(5.3) 


The expression for can be simplified by taking advantage of the 
smallness of the ratio m e /m^. If condition 

Tt «hx t Tc/m., # (5.M) 


is satisfied, formula (5*3) implies 


/20 
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( 5 . 5 ) 



In the general case, by order of magnitude 


•’£"<£ - f : ( /m.,J ( Ti /T c } fti i /) r n <t t 


( 5 . 6 ) 


Thus, when condition ( 5.*0 is satisfied, which holds nearly always 
in ionosphere conditions where the temperature of ions is lower 
than the temperature of electrons T ee << Tg i Is valid. If, in addi- 
tion, the condition 

T l «(»n l /m e ; V3 T e (57) 

is satisfied, (which is valid in the ionosphere), also holds. 

In the case when T =T, , then for example, for fully ionized plasma 

C X 

consisting of electrons and protons, we have 

Thus, local thermal equilibrium within each plasma component of the 
ionosphere is established faster than among all components of the 
system. As a result, a quasi-equilibrium state arises, which, in 
general Is characterized by three temperatures: an electron T , 

ion ^ and neutral T n temperature. This circumstance allows to de- 
rive macroscopic transfer equations with various temperatures in a 
description of physicochemical processes taking place in the iono- 
sphere of a planet. 

Elaborate thermal equilibrium of the ionosphere, and individually 
the electron and ion temperature, depends on sources and sinks of 
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thermal energy of Ionosphere particles which will be discussed In 
the next section. 

6, Equations of Heat Balance In Ionosphere 

One of the main difficulties of planetary aeronomy Is an ade- 
quate description of Inflow of heat to the upper atmosphere. By 
analogy with Earth's thermosphere, we start out with the fact that 
the energy balance of large-scale dynamic systems in the upper 
atmosphere of a terrestrial planet is essentially determined by 
absorption of short wave solar radiation by atmospheric components, 
absorption (emission) of thermal radiation from the surface of the /21 
planet and atmosphere, aeronomlc reactions, dissipation of particular 
types of waves (acoustic, gravitational and magnetohydrodynamic), 
dissipation of turbulent energy (in lower thermosphere) and Joule 
dissipation of ionosphere current heated by corpuscular flows having 
different characteristics and also dynamic processes resulting in 
a redistribution of heat from nonunlformly distributed sources. The 
general form of the intrinsic energy balance equation of a laminar 
flow of a multicomponent gaseous mixture in the upper atmosphere can 
be written as follows [12] 


*d§ 




( 6 . 1 ) 


T*i p 

Here e~F, s s 13 the thermodynamic intrinsic energy per unit volume 

of the gaseous mixture, e is the intrinsic energy per one molecule 
| ^ a 

of type s, q + is the heat flow, q^ is the inflow of energy due 
to radiation and Q d are possible local heating sources of the at- 
mosphere. The last term in the right member of equation (6.1) for 
the external mass force (1.3) is equal to and corresponds to 

Joule beating of the atmosphere by ionosphere electric currents. In 
the general case the energy inflow due to radiation consists of 
an irfrared radiation flow from the Sun and the atmosphere near the 
planet, (q IR ) and a flow of ultraviolet and X-ray radiation from the 
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Sun, leading to heating of the atmosphere during photoionl&atlon, 
photodiaooclatlon and the following aeronomic reactions. In the 
case of a two-dimensional model of the atmosphere, we have for the 
amount of short wave solar radiation absorbed per unit volume of the 
medium per unit time the expression 



(6.2) 


w: ere F Is the amount of Incident flow of solar photons on the 

yoo 

upper boundary of the atmosphere, 0 & Is the Sun’s zenith angle and 

6 (v) Is the radiation absorption cross-section of the s-th com- 

8 

ponent of the atmosphere. 

The absence of thermal equilibrium In the Ionosphere of a 
planet, In the general case, does not allow using only equation 
(6.1) written in terms of a single kinetic temperature of the gas 
during calculation of the heat balance. Separate heat balance equa- /22 
tlons muBt be used for electrons, ions and neutral particleB whose 
temperatures differ sharply under ionosphere conditions. In the 
monograph by Ivanovskiy et al [3*i], the equations of motion for the 
energy of electrons and ions are derived on the basis of kinetic 
theory. In this study we will derive these equations heuristically , 
considering the gas as a continuum. Using the definition of intrinsic 
energy of a medium e and continuity equation (3.10), we write (6.1) 
in the form 


( 6 . 3 ) 


where h =e + p /n„ is the enthalpy per one molecule of the s-th 
s s s s 

component, the coefficients y indicate the portion of absorbed 

s 

short wave solar radiation converted directly into intrinsic energy 
of the s-th component, gY s = 1, e B ^ s *^' ls the rate at which the thermal 
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energy of component s increases per unit volume because of the 
electromagnetic field. 

Assuming that viscous energy dissipation and radiant heat ex- 
change are essentially determined by the neutral component of the 
thermosphere, using the energy equation form (6.1) for the entire 
plasma, we will set up an equation which is analogous in structure, 
expressing the lav; of conservation of individual charged components 
of the ionosphere 


d 

cLt 


( S s ) =- v$ s - n s | E se 


livi 

Kit 


(6.4) 


where the last additional, term in the right member of the equation 
represents the exchange rate of energies between the Isolated s-th 
and all other components (this term drops out from the full energy 
equation (6.3), E s g 0). This term takes into account, for instance, 

the energy lost by component s (per unit time per unit volume) due 
to elastic and inelastic collisions leading to rotational, vibra- 
tional and electron excitation of neutral particles with other one- 
component media, having in the general case a temperature which is 
different from T . 


The thermal energy losses of electrons per unit volume due to 
elastic collisions with heavy and neutral particles can be determined 
using formula N N 

£j = Jj 5j 2 K (Te . (6.5) /23 


In fact, it is known that during a collision of a light particle with 
a heavy particle at rest, the portion of transferred energy may be on 
the order of magnitude of the ratio of the masses m. /m . . For example, 
during isotropic scattering, the average portion of transferred ener- 
gy is 2 m /m.,- Thus, the mean energy lost by one electron in one 
e j 
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elastic collision can be expressed in the form 


~7?rf 0 #x(T e -Tj;. 

The energy transfer rate for an electronic gas interacting with 
particles of component j is determined from the formula 





( 6 . 6 ) 


Because of the factor 2rn /nu in the right member of (6.6), this 

® J 

specific electron cooling rate (electrons undergo cooling when 
T n > T,; when T, > T , the reverse process takes place) assocla- 

S J J ® 

ted with elastic collisions is rather small (in the ionosphere 

2nu/nn ^ (2-3) *10 ). Electron energy losses occur much more rapidly 

“ J 

as a result of inelastic collisions, which were taken into account 
by introducing the Inelastic collision energy loss coefficient 
in equation (6.5) which may differ from unity by several orders of 
magnitude. In the general case, calculation of 6^ requires a know- 
ledge of the cross-sections of all inelastic collision processes 
existing in the ionosphere (vibrational, rotational and optical 
energy level excitation processes). 

When the temperature of electrons is not high (T ^ 10 J K), the 

© 

main part is played by excitation of molecular rotational energy 
levels. According to [35] > for nonpolar Ng and 0 ^ molecules, we 
have as the portion of energy <5 lost to excitation 

^<5(1*0,027^)/^ $ ao: 3,3/(i*<W6l% (6 7) 

At T > 10 K, the excitation of vibrational energy levels and at 
> ® "3 ij 

T g 'v (3*10-10 )°K, the excitation of optical energy levels already 
plays an Important part. In particular, according to [36], 
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the optical energy losaes in atomic oxygen are 


<£„ -S2,S/TT. V * 


( 6 . 8 ) 


and owing to excitation of transitions between fine structure levels, 

they become noticeable at temperatures which are characteristic for 

ionospheres of planets (T 'v 10 Jo K). 

© 

For a mixture of gases, S Is determined according to formula /2*l 





( 6 . 9 ) 


where 3-„ and 6 are the collision frequencies and the portion of 
transferred energy for a gas of type s. Mention should be made of 
the fact that the method of taking Into account the inelastic col- 
lision energy loss during collisions of electrons with molecules 
in terms of the coefficient 6 Is approximate even in the case of 
considerably high energy transfer rates between rotational and vibra 
tional degrees of freedom on one hand and translatory degrees of 
freedom on the other hand. A detailed calculation is necessary when 
the relationship is weak [37] . 

In the case of interaction of heavy particles i and n we have 

[ 18 ] 


from which follows the energy transfer rate for an ion colliding with 
neutral particles n 


_ - m. m. 



( 6 . 10 ) 
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Let us now write dov. r n the energy balance for an electron gas in 
the case of three-temperature plasma. For electrons we have 
£ B 3/2 KT q and h Q =5/2KT , so that equation (6.4) after formula 
(6.2) has been taken into account assumes the form 



( 6 . 11 ) 


where Q dg = 2-lcT 12 <f> e [eV/cm 2 • see ] (4> e represents a flow of fast 
electrons with energy E) is the Inflow of heat to the thermalized 
electron gas of the ionosphere from high-energy particles which 
fell out into the atmosphere [38] , L are the heat losses of elec- 

C 

trons during rotational excitation of nonpolar molecules, during 
vibrational excitation of molecules, and during excitation of tran- 
sitions between levels of the fine structure of atomic oxygen [39]. 
Since Joule heat is essentially released in electron gas, to empha- 
size this fact, it is convenient to represent the term -eJ e -S' in 
the form f c £ _ where J is the current flowing 

in the ionosphere as a result of the action of an electric field 
(for example auroral current or an equatorial electro jet in the 
Earth's ionosphere). 

An additional source of heat for the ionosphere of a planet, /25 
especially in the case of nonmagnetic planets, may be hot solar wind 
plasma, the interaction of which with the planet constituting an 
obstacle leads to generation of hydromagnetic and acoustic energies. 
Another heating source in the topside ionosphere of a planet is the 
turbulent heat conductivity of hot solar wind plasma. We will pre- 
sent some estimates : 

The energy transfer rate E e , n for an electron colliding with 
neutral particles, for the principal neutral components of the 
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Earth’s ionosphere has the following values: 


Ev.-I.77 

ro ^ /\ ue - TV eV/ cur* sec E^ «*l . 21 • icfto, «.<> t /I+3 . 6 • IO l T« w / 
T.^Te-T^ i E e ,a =3.74 .itfTfcrf.V'V WnJi E v ,=9.63.I0‘'V e 
/I-I.35 10 'VWTWk.J { ^*-8.46 lO m %\T^(T^T K )[2l,2z\, 


The rate at which an electron gas loses energy which is trans- 
ferred to a mixture of positive ions, for example ot He and H + is 
defined in [40] as 


■[ 40 ] »' 

A E, 

Jife 


_&1Q 


°i 


eV/cn,3 * sec (6.12) 


According to [40] cooling of electron gas by rotational exci- 
tation of N g and 0 2 molecules can be written in the form 


L^^.^IO L f ^6 1 9 • 10 /7 C 7c 


t/i 


(%-V. 


(6.13) 


At T > ^1500°K, vibrational excitation [4l] is more effective for 
N 2 and 0 2 


eV/crn 3 .sec 


M'-to exp[$ aoooi^ 
y = f •/. iO Hk 10 ’ i 0 ~ 3 ( T «- 18 0 0) j ^ / £ 1 V 

f «3500+-/,£33flW<>0c>)-2, OS6'f0^(Tt-10Qo)(T e -‘toao); 

\! %h ^B-io- v n^^^\^ooo^X 

9 '« 3,3 0 2 - ■ 10 3 t h, 3 $ '1 o' 3 L fi [ k',56' • 10 ' " H T e - 2 4 0 0 ) ] t 


(6.14) 
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Losses by excitation of the fine structure of the 0 atom ground 
level and the first excited metastable J 'D level of the same atom can 
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be written respectively in the form 


\! Q0 n~Z t ^iO m1l tt a lt 9 (i-fr10 m \)(T 6 - T*,)T*\ eV/cm 3 -sec 

~(q$i 3+0, f&3' to %)<*£- J •{tytSi+qwt'io 

*1 


Calculations of y^ during high and low solar activity were car- /26 
ried out in [43-453. According to [453, we have for y^ 


_./ f eV for 12.04 600«h. 

^[•/O V 3 ° eV for 7£6O0*«. 


( 6 . 16 ) 


We will write the intrinsic energy transfer equation for the 
ion component of the atmosphere ignoring viscous dissipation of ions 
as follows 



(6.17) 


where 



is the intrinsic enero-v ner unit volume of the ion comno- 
nent cf the ionosphere. Heating of ions takes place during Coulomb 
Interaction with electrons, in this case, the rate at which heat 
arrives in the ion gas is equal to the cooling rate of elec- 

trons in equation (6.11). tfhe ions undergo cooling during elastic 
collisions with heavy particles. In this case the heat transfer 
rate is either determined by the difference between the temperatures 
of the ion and neutral gases or by the possible difference in the 
directed velocities and (frictional heating). Effective cooling 
of ions in the "parental" gas takes place during resonance transfer 
of the charge; In this case, the expression for is changed by 
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the factor ( T n + ^ ion ) > where 0.3 ^m^ n.ij for most atmospheric gases 
[^6], For the Earth’s ionosphere, the heat losses consist of elas- 
tic collision losses of 0* He^ H ions (taking into account resonance 
recharging) [21,22] 



( 6 . 18 ) 


The inflow of heat to the ion gas from photoionization is small 
and it is not taken into account in equation (6.17). Usually the 
Ion temperature is determined from the equilibrium form of this 
equation, while the effect of the thermal conductivity of ions com- 
pared with the thermal conductivity of electrons (in view of the 
smaller thermal conductivity rates of ions) can also be Ignored, 
since X e /Xj ^ /m e /m 1 >> 1 [ 18 ], 

Kinetic theory of multicomponent gaseous mixtures of single- /27 
atom gases of moderate density in any approximation obtained by the 
Chapman-Enskog method leads to the following expressions for heat 
flows transferred by electrons and ions [18,293 


* «' A **■'£«* (6.19) 

A vl^ “ 2 Ji , 

f a i*V * -» 

Aie 2 7 Turn Je ~ / 

in 4 (6.20) 
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A A A A 

where X , X. and X , X are the thermal conductivity and mutual 

thermal conductivity tensors for electrons and ions,/f e , ^ and 
<^ > ei J A'le are thermoe ^ ec ^ r ^ c force tensors. Usually the contri- 
bution of cross-effects to the total rate of the process is smaller 
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by one order of magnitude in comparison with the direct effect [12], 
Moreover, at a moderate degree of ionization of Ionosphere plasma, 
when collisions with neutral particles play a crucial part, the 

A A 

reciprocal heat conductivity tensors X e ^ and X ie and thermoelectric 
force tensors c/f and j/e are Insignificant [29], The second sim- 
plification is related to the circumstance that the heat flow is 
mainly determined by the thermal conductivity XAt, since on the 
strength of the quaslllnearity of plasma, the ion and electron dif- 
fusion flows have the same order of magnitude and thermal conductivity, 
not coupled by quasilinearity conditions proceeds much faster 
( times faster). 

Thus, when S-0, formula (6.19) for a heat flow transferred by 
electrons has the form [18] 


t Vc--^e7Ta ^ 


* _ <?, *# K a b c Te 


( 6 . 21 ) 


where X is the total electron thermal conductivity of the gaseous 
e 

mixture and is the total mean collison frequency of electrons with 
all heavy particles. When 5^0, formulas (6.21) are valid for a heat 
flow which is parallel to the magnetic field. There formulas take 
on a different form for the perpendicular component of the heat flow. 
[18]. 


In the topside ionosphere of planets with a magnetic field, 
heat transfer in electron gas occurs only along the lines of force 
of the magnetic field (the ratio Xj,/X jj is an extremely small quan- 
tity, and essentially, it depends on electron-ion collisions. In 
this case, thermal conductivity is the same as that for a completely 
ionized gas [30, ^7] 


/ 




Soffit (2 
e*firT e €hA*i 


( 6 . 22 ) 
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Since the atmospheres of planets represent Incompletely Ionized /28 


plasma, especially at low altitudes, where considerable collisions 
occur between electrons and neutral particles, in the general case 
expression (6.22) is not valid and formula (6.21) must be used which, 
taking into account (6.22), can be represented in the form 


X 

Ae 


1 

Afll $ A«k ; 


(6.23) 


where is the thermal conductivity of electrons in a neutral gas 
composed of various components (collisions of charged particles are 
ignored). This quantity is given by 






-l/Sw/ffi: V^i 
3 W'VflCmJ Q ck 


(6. 24) 


The thermal conductivity coefficient of tvDe i ions is 





(6.25) 


Using the full energy equation and heat balance equation for 
an electron and ion gas, an energy equation can be derived for a 
neutral atmosphere with heat sources (sinks) because of the inter- 
action with the plasma component of the planet's ionosphere. These 
equations permit to investigate thermal conditions in a neutral 
atmosphere and ionosphere plasma, while taking into account the 
mutual effect of components of the system. 
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